Light clusters and the pasta phase 
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The effects of including light clusters in nuclear matter at low densities are investigated within 
four different parametrizations of relativistic models at finite temperature. Both homogeneous and 
inhomogeneous matter (pasta phase) are described for neutral nuclear matter with fixed proton 
fractions. We discuss the effect of the density dependence of the symmetry energy, the temperature 
and the proton fraction on the non-homogeneous matter forming the inner crust of proto-neutron 
stars. It is shown that the number of nucleons in the clusters, the cluster proton fraction and the 
sizes of the Wigner Seitz cell and of the cluster are very sensitive to the density dependence of the 
symmetry energy. 
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I. INTRODUCTION 

The formation of light clusters in nuclear matter at low 
densities and its influence in the appearance and compo- 
sition of pasta structures has been frequently discussed 
in the literature |l|-KL(j|. 

Below saturation density homogeneous nuclear matter 
can become unstable against phase separation and sev- 
eral types of complex structures can be formed as a re- 
sult of the competition between the strong and the elec- 
tromagnetic interactions. This pasta phase |lll - ll5i is 
found at densities of the order of 0.001 - 0.1 fm~^ [l^ in 
neutral nuclear matter, formed by protons, neutrons and 
electrons, and in a smaller density range in /3-cquilibrium 
stellar matter [l3,[3- ^t very low densities, up to 0.001 
times the saturation density, and moderate temperatures, 
the few body correlations are still important and the 
system minimizes its energy by forming light nuclei like 
deuterons, tritons, helions and/or a particles. The ap- 
pearance of these light clusters can modify the behavior 
of the neutrinos in the expanding matter resulting from 
a supernovae core collapse and affect the cooling process 
of the protoneutron star. 

In a previous paper [l9| we have studied the influence 
of the a particles both in homogeneous matter and in the 
onset and structure of the pasta phase within the rela- 
tivistic mean field approximation. We have considered 
both free a particles and a particles interacting through 
an a — a; meson coupling. This repulsive interaction is 
essential to avoid an overprediction of a particles above 
p ~ 0.001 fm~'^ and it is also the mechanism responsi- 
ble for the dissolution of the a particle clusters in the 
approach considered. 

Model dependences are more important for p > 0.001 
fm""^ when the a particle fraction differences between 
models may be as large as one order of magnitude or 
even larger. The effect of the temperature is to shift the 



maximum of the a particle distribution and the density of 
cluster dissolution to larger densities. At the same time, 
the maximum values of the distributions themselves were 
shown to decrease with the increase of the temperature. 
The maximum values of the a particle distributions also 
decrease when the proton fraction decreases. However, 
the proton fraction has no effect neither on the density 
localization of the maximum, nor on the density of dis- 
solution of the clusters. 

It was also shown that in the pasta phase formed in 
asymmetric nuclear matter the a particle fraction in- 
creases with temperature. This is an interesting effect re- 
lated to the proton fraction in the background gas which 
increases with temperature for asymmetric matter. It is 
important to test the above mentioned behaviors when 
other light clusters are also included in the system. 

In the present paper we extend our previous work [l9| 
by considering also deuterons, tritons and helions. We 
study the distribution of these light clusters as a func- 
tion of the baryonic density for several temperatures and 
proton fractions, and investigate the effect of the light 
clusters on the onset, type of structures and size of the 
heavy clusters {A > 4) of the pasta phase. We compute 
the mass and charge content of the droplets formed in 
the droplet regime of the pasta phase and compare with 
other approaches. 

Our approach to the description of the pasta phase is 
within what is known as the single nucleus approxima- 
tion. The same approach was used in 1, 2], however, only 
alpha particles were included as light clusters. In [3] a 
statistical model consisting of an ensemble of nuclei and 
interacting nucleons in nuclear statistical equilibrium was 
proposed to describe supernova matter. There, it was 
shown that the presence of light clusters besides alpha 
particles are of particular importance at low densities. A 
similar conclusion was drawn in [8| where a phenomeno- 
logical statistical model consisting of free nucleons de- 
scribed within a mean-field approximation and a loosely 
interacting cluster gas was formulated to describe mat- 
ter in supernova explosions and proto-neutron stars. In 
all these works an excluding volume concept was used: 
in the first two works this was done with respect to the 



alpha particles and in the other two with respect to all 
nuclei or clusters. 

The validity of the two approximations, excluded vol- 
ume and single nucleus has already been discussed. Re- 
cently, the excluded volume approach was compared to 
two quantum many-body models in [10[ and it was shown 
that this approach is a bit crude at temperatures of the 
order of 5 MeV, although this occurs at densities where 
the composition of matter is dominated by heavy ions. In 
j2(j | it was discussed that the EOS is not much affected by 
the single nucleus approximation, which describes matter 
composition in an average way, although, a correct distri- 
bution of nuclei may be important to describe correctly 
the supernova-dynamics. 

The approach we consider in the present work has the 
drawbacks of a single nucleus approximation and does 
not include shell effects contrary to statistical models 
[l^. Therefore, we will restrict ourselves to tempera- 
tures above which shell effects are no longer important. 
However, we would like to point out that at densities 
close to the crust-core transition statistical models have 
difficulties in describing properly the medium effects on 
the nuclei and exotic structures such as the pasta phases 
contrary to the present approach. 

As referred before, we will avoid the excluded volume 
approximation and will include Pauli blocking and self- 
energy effects in a phenomenological approach by includ- 
ing a meson-cluster interaction within a relativistic mean- 
field formalism. A quantum statistical approach allows 
the calculation of the medium infiuence on the quasi ener- 
gies of clusters @,[a,lMl- I^i particular, recently analytical 
fits to the quasi-energy shifts of light nuclei were pub- 
lished Q . These results allow the determination of a bet- 
ter parametrization of the cluster-meson interaction than 
the one used in the present work. However, we should 
point out that when the medium effects are more im- 
portant the contribution of light clusters is already very 
small. Moreover, we have shown in [l^l that the fraction 
of alpha particles obtained with the prescription used in 
the present work agrees with the results both of the virial 
equation proposed in [J] and the results of [6| within a 
generalized mean-field model. 

A relativistic mean-field approach is a phenomeno- 
logical theory where the meson-nucleon and the meson- 
meson interactions mimic the different contributions of a 
quantum many-body formalism such as the Pauli block- 
ing and self-energies and it is difficult to separate the 
different contributions which determine the behaviour of 
the system. The same may be said with respect to the 
meson-cluster interactions introduced in [y, ll9| . In a re- 
cent study [21| the effect of the meson-cluster couplings 
in the equation of state (EOS) of homogeneous nuclear 
matter with clusters was investigated within a zero tem- 
perature relativistic mean-field approach. It was shown 
that the tr-cluster and w-cluster couplings determine the 
behavior of the clusters, namely their fraction and disso- 
lution density. 

In this work we use the relativistic mean field ap- 



proximation and we consider four different parametriza- 
tions. We have chosen the NL3 [HI , NL3wp ^ , FSUgoW 
[2JI and lU-FSU [23 parametrizations of the non-linear 
Walecka model (NLWM) [2g|, which allow us to discuss 
the role of the density dependence of the symmetry en- 
ergy on the properties of the non-homogeneous nuclear 
EOS. 

The paper is organized as follows: in section II we 
briefly review the formalism underlying the homogeneous 
neutral npe matter with the inclusion of the light clus- 
ters. In section HI the coexisting-phases method used to 
obtain the pasta phase is briefly reviewed. In section IV 
our results are displayed and commented and in section 
V the final conclusions are drawn. 



II. THE FORMALISM 

We consider a system of protons and neutrons with 
mass M interacting with and through an isoscalar-scalar 
field (j> with mass nis, an isoscalar- vector field V^ with 
mass TUy, an isovector- vector field b^ with mass nip. We 
also include tritons ('^H, represented by t), helions ("^He, 
represented hy h), a particles and deuterons (d). A sys- 
tem of electrons with mass nie that makes matter neutral 
is also included. 

The Lagrangian density reads: 

+ C^p + Ce+CA, (1) 

where the Lagrangian density Cj is 

Cj^^,j[jpiD^^M*]^,, (2) 

the a particles and the deuterons are described as in @ 
with Ca and Cd given, respectively, by 



and 



^d = ^(*-Dd</'d - iDd'PdT^'^Ddiii'du - iDdu4'd,i) 



hi^* nr'^^ 



^ 2^d Md(j)d^i, 



(4) 



with 

-J — -- ^vj ' 2 

3 =P,n,t,h,a,d 

M* = M - g^cj), j =^p,n 

M* = Mt ^MI - Bt, 

M* = A4 = 371/ - B„, 

M* = Md^2M~Bd, 



d^ - gv,V^ -^T-h^--{l+ T3) A^ (5) 



(6) 
(7) 
(8) 
(9) 
(10) 



with the binding energies given by Bt = 8.482 MeV, 
Bh = 7.718 MeV, B^ = 28.296 MeV and B^ = 2.224 
MeV and g^j ~ Ajg^ and gpj ~ \Zj — Nj\gp, where Aj 
is the mass number, Zj, the proton number and Nj, the 
neutron number. Notice that in our model the cluster 
masses are fixed to constant values, they do not depend 
upon the temperature and the density. However, due to 
the cluster-meson interaction, we find (see Sec. IV) that 
the dissolution density of each cluster increases with the 
temperature. 

The electron Lagrangian density is given by 



and the remaining terms in eq. ([T|) are 



(11) 



Ca = 

r — 

^P = 

C — 

Ca = 



1 



d. 



- ml<i? - \k^^ - ^M' 



\ l-Ui^,^^'^ + ^nlVpV^ + ^^igliy^v^f 



"p^p 



t^.glglV^y^\y^. ■ b^ 



1 



F„„F^'' 



where £1^^ = d^V^-d^V^, 'Bpy = di_ihy-di,h^-T p{hp x 
bj,) andF^^ — dp_A^—d,yAp^. The parameters of the mod- 
els are: the nuclcon mass M = 939 MeV, the coupling pa- 
rameters gs^ gv, gp of the mesons to the nucleons, the self 
interacting k, A and S, constant couplings and the to — p 
coupling A„. In the above Lagrangian density x is the 
isospin operator. When the NL3 parametrization is used, 
Ay is set equal to zero. When the NLSwp parametriza- 
tion is chosen, we follow the prescription of [23|, where 
the starting point was the NL3 parametrization and the 
gp coupling was adjusted for each value of the coupling 
Au in such a way that for kp = 1.15 fm~^ (not the sat- 
uration point) the symmetry energy is 25.68 MeV. In 
particular, in this work we set A„ to a moderately large 
value (A„ — 0.03, sec Table [l|. The comparison of NL3 
and NL3ajp results is meant to show the effect of changing 
the isovcctor part of the Lagrangian only. 

The FSUgoU [2J| parametrization was chosen because 
it has the advantage of reproducing some of the results 
[27| obtained with more sophisticated density dependent 
hadronic models [2^ without the need of rearrangement 
terms. FSUgoM combines the inclusion of the vector self- 
interaction term present in Ref. |29l |30[ that is respon- 
sible for explaining some observed properties at nuclear 
density and the isoscalar- isovector coupling present in 
[23| , capable of improving the density dependence of the 
symmetry energy. An interesting comparison of the re- 
sults for the neutron star mass-radius relation involving 
NL3 and FSUgoM is given in [31|. We also consider the 
lU-FSU parametrization which keeps the main properties 
of FSU but was readjusted in order to allow for neutron 
star masses up to about 2.0 Mq ^^. The parameter 



sets for the NL3, NL3wp, FSUgom and lU-FSU models 
are shown in Table HI Their corresponding bulk nuclear 
matter properties are given in Table |TT1 In order to clarify 
the discussion, we show in Fig. [TJ for the models consid- 
ered, the symmetry energy at densities below 0.3 fm^'^ 
and the surface tension at T = 5 MeV as a function of 
the proton fraction. The surface tension was determined 
according to the parametrization given in eq. (j33p and 
the Appendix. With the present choice of models we 
will be able to discuss the implications of the symme- 
try energy on the pasta phase. One expects, generally 
speaking, two types of effects: a smaller L corresponds 
to a larger surface tension for asymmetric matter [see 
Fig. [Hd) and Table HI] : a larger Sgym. leads to a more 
isospin-symmetric liquid phase. We will see in section 
IV that in models with a larger surface tension the pasta 
phase sets in at higher densities and the drip of particles 
is unfavored, giving rise to a lower density background 
gas. 



TABLE I: Parameter sets for the models used in this work. 
The masses of the mesons are in MeV and the other quantities 
are adimensional {k is given in nucleon mass units). 





NL3 


NL3a;p 


FSUcoid 


lU-FSU 




[2^ 


[23] 


[24J 


[25] 


rus 


508.194 


508.194 


491.500 


491.500 


rriv 


782.501 


782.501 


782.500 


782.500 


nip 


763.000 


763.000 


763.000 


763.000 


9s 


10.217 


10.217 


10.592 


9.971 


9v 


12.868 


12.868 


14.302 


13.032 


9p 


8.948 


11.2766 


11.767 


13.590 


K 


4.384 


4.384 


1.7976 


3.5695 


X 


-173.31 


-173.31 


299.11 


2.926 


c 


0.00 


0.00 


0.06 


0.03 


A„ 


0.00 


0.03 


0.03 


0.046 



TABLE II: Nuclear matter properties at the saturation den- 
sity and zero temperature: binding energy per nucleon B/A, 
density po, eflective nucleon mass M* , incompressibility K, 
symmetry energy Saym. , and slope L of the symmetry energy 

Osym. \P) . 





NL3 


NL3ajp FSUgoM 


lU-FSU 




[22] 


[23] 


m 


[25] 


B/A (MeV) 


16.3 


16.3 


16.302 


16.4 


po (fm-^) 


0.148 


0.148 


0.148 


0.155 


M*/M 


0.60 


0.60 


0.62 


0.62 


K (MeV) 


272 


272 


227.9 


231.2 


£aym. (MeV) 


37.4 


31.66 


32.54 


31.3 


L (MeV) 


118.32 


55.23 


60.39 


47.2 
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FIG. 1: Comparison of the a) symmetry energy, b) surface tension at T = 5 MeV calculated with the models NL3, NL3a;p, 
FSUgow and lUFSU. 



After this discussion that motivates our choice of 
parametrizations, we go back to the description of the 
method. From de Euler-Lagrange formahsm we ob- 
tain coupled differential equations for the scalar, vector, 
isovector-scalar, nucleon and cluster fields. In the static 
case there are no currents and the spatial vector com- 
ponents are zero. Moreover, we neglect, as usual, the 
Coulomb interaction in the case of homogeneous matter. 
In the calculation of the pasta phase, its effect on the 
protons will be included. In |16| a complete description 
of the mean-field and Thomas-Fermi approximations ap- 
plied to different parametrizations of the NLWM is given 
and we do not repeat it here. The equations of motion 
for the fields are obtained and solved self-consistcntly and 
they can be read off [1^, Il3 ■ The above mentioned equa- 
tions of motion depend on the the equilibrium densities 



P = Pp + Pn + 4pa + 2pd + 2,pt + iph, 



P3 = Pp - Pn ~ Pt + Ph, 



Ps 



Ps„ 



(12) 



(13) 



(14) 



The quantities pa and pd are discussed next, the 
fermionic densities are 



-^ I P^dp{f,+ - f,^), i=p,n,t,h 



and the corresponding scalar densities are 



Ps. = ^ I P^dp'-^{f,+ + fy 



m; 



(15) 



(16) 



The distribution functions are given by 
1 



/.± = 



l + exp[(e*(p)T^.O/r] 



, i ~ p,n,t,h (17) 



where e* = 
tentials are 



p^ + M*", and the effective chemical po- 



T /- dpi I 

9viVo - -Y ^3i bo, 



where 



Pt=Pp + '2pn, Ph='2pp+ Pn, 



(18) 



(19) 



and T^i = ±1 is the isospin projection for the protons 
(helions) and neutrons (tritons) respectively. 

In the present work the a particles and the deuterons 
are included as bosons and their chemical potentials are 
obtained from the proton and neutron chemical poten- 
tials by imposing the chemical equilibrium, as in [6| : 

Pa^2{pp + pn), Pd^Pp+Pn- (20) 

Their effective chemical potentials read 

^i = Pj - 9jVo, J = a, d. (21) 

The density of thermal a particles and deuterons arc 



Pj 



^ j P'dp{fj+ 



fj-) 



with the boson distribution function given by 

1 



/.± = 



-l-Hexp[(e,T^,)/T]^ 



(22) 



(23) 



where e, = 



pz 



Mf. We should point out that at low 



enough temperatures the as and deuterons contribute 
with two terms: a condensed fraction and a thermal 
contribution. Both contributions can be included ex- 
plicitly in the present formalism. We only indicate the 
thermal contribution because we have limited our discus- 
sion to temperatures above Bose condensation of as or 
deuterons. 



For the free electrons, their density and distribution 
functions are the same as for the other fermions, where 
/ie is the electron chemical potential and Ce = i/p^ + mf . 
We always consider neutral matter and therefore the elec- 
tron density is equal to the total charge density of the 
charged particles. 

In the description of the equation of state of a system, 
the required quantities arc the baryonic density p, energy 
density £, entropy density <S, pressure P and free energy 
density T = E — TS, and their expressions are explicitly 
given in [11 [11. 



III. NUCLEAR PASTA 

In this section we describe briefly the coexisting-phases 
(CP) method to study the non-homogeneous phase of nu- 
clear matter with a fixed proton fraction. The basic idea 
is that matter can be organized into separated regions 
of higher and lower density, and the geometry of these 
regions is assumed to be very simple: a lattice of spheri- 
cal droplets (bubbles), a plan of cylindrical rods (tubes) 
or an alternating sequence of slabs. The interface be- 
tween regions is sharp, and it is taken into account by a 
surface term and a Coulomb one in the energy density. 
In the spirit of this approach, a single geometry will be 
the physical one for some given conditions (temperature, 
density and proton fraction). 

As in [15|, [la] , for a given total density p and proton 
fraction, now defined including the protons present in- 
side the light clusters, the pasta structures are built with 
different geometrical forms in a background nuclcon gas. 
This is achieved by calculating from the Gibbs conditions 
the density and the proton fraction of the pasta and of 
the background gas, so that we have to solve simultane- 
ously the following equations together with eqs. dH]) and 
(EOj): 



pi = p 



II 



p.j — fij^ , i = p,n,t,h,a,d, 

I pi + (1 - f)p" = yp9 P: 



(24) 



(25) 



(26) 



where I and II label the higher and the lower density 
phase respectively, and / is the volume fraction of phase 



^ P'-P" 



(27) 



The total baryonic density p is given by eq. (IT2]) . Y^g is 
the global proton fraction given by 



V - ^ 



(28) 



and Pc = Pp -^ '2pa + Pd + Pt + '^Ph stands for the charge 
density. 



The density of electrons is uniform and taken as Pe = 
Ypg p. The total pressure is given by the sum of the 
nuclear and electron partial pressures P = Pnuci + Pe, 
where the nuclear contribution Pnuci includes the kinetic 
contribution of each type of particle (nucleons and light 
clusters) plus the meson contribution. The total energy 
density of the system is given by 



£^f£' + il- f)£" + £e + Esurf + £, 



Coul 



(29) 



By minimizing the sum Ssurf + £coui with respect to the 
size of the droplet/bubble, rod/tube or slab one gets [15 1 
Esurf = 2£coid, and 



2F 
Seoul = ji7j(eV$) 



1/3 



[^D{pi-pi')]- 



(30) 



where F is the volume fraction of the inner part (F — f 
for droplets, rods, slabs and F = 1 — f for bubbles and 
tubes), a is the surface tension, D is the dimension of 
the system, and $ is a coefficient depending upon F and 
D 0,[ll. 

Each structure is considered to be in the center of a 
charge neutral Wigner-Seitz cell constituted by neutrons, 
protons and electrons [2]. The Wigner-Seitz cell is a 
sphere/cylinder/slab whose volume is the same as the 
unit BCC cell. In [2| the internal structures are associ- 
ated with heavy nuclei. Hence, the radius of the droplet 
(rod, slab), Rd, and of the Wigner-Seitz cell, Rw are 
respectively given by: 



R 



D 



aD 



4^e2(^/_^//)2$^ 



1/3 



R 



Rt 



w 



F^/D' 



(31) 



(32) 



Concerning the surface energy, the authors of [32| have 
shown how the surface energy affects the appearance of 
non-spherical pasta structures. Also the authors of [l^ 
state that the appearance of the pasta phase essentially 
depends on the value of the surface tension. We have 
fixed the surface tension at different values and confirmed 
their claim in |16l . llSj . where the surface tension was 
parametrized in terms of the proton fraction according 
to the functional proposed in |33j . obtained by fitting 
Thomas-Fermi and Hartree-Fock numerical values with 
a Skyrme force. In [l9[ a new prescription was used and 
the surface energy was fitted to the results obtained from 
a relativistic Thomas-Fermi calculation and its depen- 
dence upon the temperature was taken into account. A 
mathematical formula for a that gives accurate results 
for a broad range of neutron excess and for temperatures 
up to 10 McV is 

a{x, T) = a{x) [l - a{T)xT - b{T)T^ - c(T)x^T] , 

(33) 
where a'{x) is the surface tension at T = and x = 6^ 
stands for the squared relative neutron excess: 



Pp 



1 - 2Yr„ 



In Tabic Hm the a parameters fitted to the Thomas-Fermi 
approximation results up to T = 10 MeV are given. For 
T ~ IQ MeV the parametrization is good if Ypg > 0.2. 
For lower temperatures, the range of validity extends to 
lower Ypg values. 



IV. RESULTS AND DISCUSSIONS 

In this section we present the results of this work. In 
Fig. [2] the free energy per particle is shown for four differ- 
ent constructions of the EOS: homogeneous matter (H), 
homogeneous matter with clusters (HC), mixed homo- 
geneous matter and pasta phase (HP) and finally HP 
with clusters (HPC). Notice that with HP we mean that 
both non-homogeneous matter in five different shapes 
(droplets, rods, slabs, tubes and bubbles) and homoge- 
neous matter are computed and the phase with the lowest 
free energy per particle is the physical one. 

The free energy is clearly lowered by both the inclusion 
of pasta and clusters. However, the pasta phase affects 
the results in a large density range and leads to a decrease 
of a few MeV [Fig. 2a)], whereas the clusters are visible 
only at very low density and lead to a decrease of ^ 1 
MeV [Fig. 2b)]. Their effect is even smaller in the pasta 
phase range. It is also seen, as already discussed in Q, 
that the pasta phase sets on at a larger density when 
light clusters are included. 

The inclusion of the light clusters in the EOS improves 
the CP method. In fact the way the surface energy en- 
ters in the CP calculation doesn't allow the appearance 
of these light nuclei but, as we have seen, their presence 
is important just below the onset of the pasta phase. We 
will later verify whether the inclusion of the light clusters 
explicitly does not correspond to a double description of 
the small clusters in the pasta phase. In fact, if we in- 
clude light clusters and obtain droplets formed by four or 
less nucleons the same physical object is being described 
simultaneously in two distinct ways. 

In the more complete case (i.e., HPC), we can see how 
the fraction of nucleons belonging to different structures 
changes with the total baryonic density. Nucleons can be 
free (not clusterized) , belong to a light cluster {A < 4) 
or to a heavy one. In our formalism, the role of heavy 
clusters is played by the high density part of pasta struc- 
tures. In the regimes of temperature and proton fraction 
where a non- homogeneous phase appears, the following 
sequence is generally found: homogeneous matter at very 
low baryonic density, pasta phase at intermediate densi- 
ties and then again homogeneous matter at higher den- 
sities. 

In the non-homogeneous matter range, the fraction of 
a generic constituent has a contribution from the high 
density (phase I) and one from the low density phase 
(phase II): 



where i = p,n, a, h, t, d and / is the volume fraction of 
the denser phase defined in (P?]). 

Therefore, we classify the fraction of nucleons as fol- 
lows: 



^71, free 
flight 
^ heavy 



Yn 





Yh+Yt + Yd 



(35) 



in homogeneous matter at low density. 



^p^free 

^ light 
^ heavy 



II 



= Y,V{l-f)^ 



{Y^ 
J p 



Y^' + Y/' + Yi'){l-f)^^ 



in the non-homogeneous phase, and finally 



^ p. free 
^n,free 
^light 

^ heavy 



(36) 



(37) 



Y, = Y/f^ + Y/'{l-f)^ 



(34) 



in the high density homogeneous phase. 

In Fig. 13] the fractions of nucleons are shown for the 
NL3 and NL3a;p parametrizations at T = 5 MeV, and 
Ypg = 0.2. For both parametrizations, matter in the low 
density regime is formed by not-clusterized neutrons and 
protons and a light-cluster fraction. The magnitude of 
this fraction depends on the temperature. At p ^ 0.001 
fm""^ the pasta phase sets on. The background low den- 
sity gas is constituted mainly by neutrons with a small 
fraction of protons and light clusters. For densities above 
0.08 fm^'^ there is a transition to dense homogeneous 
matter, interpreted in the present approach as an infi- 
nite cluster. 

The differences between the results obtained for NL3 
and NLSwp can be interpreted as a consequence of the 
effect of the density dependence of the symmetry energy. 
Both NL3 and NLSwp have the same isoscalar properties. 
However they differ on the isovector channel, namely the 
symmetry energy and its slope at saturation are different, 
NLSwp having a smaller symmetry energy and smaller 
slope at saturation. Below p = 0.1 ini~^ the symme- 
try energy is larger for NL3w/9 while the opposite occurs 
above that density [see Fig. [Hd)]. It is seen that both 
models behave in the same way below the onset of the 
pasta phase: the proton, neutron and light cluster frac- 
tions are practically equal. However there are noticeable 
differences on the pasta phase. The most important ones 
are the onset density of the pasta phase and the fraction 
of nucleons in the cluster. A smaller symmetry energy 
slope shifts the onset of the pasta phase of NLSwp to 
larger densities because it gives rise to a larger surface 
energy that hinders the formation of pasta structures. A 
larger surface energy also gives rise to a smaller fraction 
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FIG. 2: (Color online) Free energy per particle as a function of the baryonic density for NL3 at T = 5 MeV and Y^g = 0.2. 
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H and HC or HP and HPC are indistinguishable; b) at very low densities. 
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FIG. 3: (Color online) Fractions of nucleons as a function of the density, for the NL3 and NL3(^p parametrizations at T = 5 



of neutrons outside the cluster because it is more difficult 
for neutrons to drip out. Since neutrons have an impor- 
tant role on the cooling of the crust, the fraction of free 
neutrons on the pasta phase range will certainly affect 
the cooling and transport properties of the crust. 

We plot in Fig. |3] the fractions of nucleons predicted 
by the FSUgoW and lU-FSU parametrizations compared 
to the NL3 results. There are some similarities between 
these results and those obtained with NLSoip. The on- 
set of the pasta phase density occurs for the lU-FSU 
parametrization at values slightly higher than in FSUgoU 
and much higher than in NL3, due to its smaller L 
and, therefore, larger surface tension. At low densities, 
the main differences among the three models occur at 
p ~ 0.001 fm~^ with different proton and light cluster 
fractions. These differences are mainly due to differences 
in the density dependence of the isoscalar channel of the 
EOS. The FSUgoM parametrization has the largest gs 



coupling, so, since nucleons have a smaller effective mass, 
and a larger binding energy per particle, the formation 
of light clusters is not so favored as in a model where the 
nucleon effective mass is larger. However, a small frac- 
tion of light clusters is still observed in FSUgoW , because 
the onset of pasta phase occurs at densities larger than 
in NL3. 

The effect of the temperature on the nucleon ratios is 
examined in Fig. [5] where the results for T = 3.5, 5, and 
6.5 MeV are shown for the lU-FSU parametrization with 
Ypg = 0.2. The main features of increasing the tempera- 
ture can be summarized in three points: a) the onset of 
the pasta phase (of the core) is shifted to higher (smaller) 
densities, as already discussed in jSJ]; b) the low density 
gas of the pasta phase has a larger number of particles; 
c) the light clusters contribution is smaller at densities 
below the pasta phase onset and larger on the low den- 
sity background gas of the pasta phase. These behaviors 
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all occur because when the temperature increases the in- 
stability region decreases and there is a larger number of 
nucleons that drip out of the dense clusters. 

The effect of the global proton fraction on the pasta 
phase is clearly seen in Fig. [SJ where the different con- 
stituent fractions at densities below 0.1 fm"'' are plotted 
for Ypg = 0.2, 0.3, and 0.5 for the IU-FSU parametriza- 
tion at T = 5 MeV. The onset of the pasta phase occurs 
at larger densities for the more symmetric matter. This 
is possibly a trend due to the method used for the pasta 
phase clusters with a zero thickness surface. In fact, in 
the pasta calculation within the Thomas Fermi approach, 
where the surface is described in a self-consistent way, 
the opposite occurs [SJI- However, if we also consider 
the light clusters, the onset of the clusterization starts 
at lower densities for the larger proton fractions. In the 
pasta phase, the larger proton fractions give rise to clus- 
ters with a larger number of nucleons immersed in a gas 
of nucleons and light clusters with a lower density. 



These results allow us to make some comments on the 
possible consequences of the evolution of a protoneutron 
star. After the supernova explosion neutrinos are trapped 
inside the star and we can consider that the lepton frac- 
tion is approximately constant, taking a maximum value 
of 0.4 and decreasing as neutrinos leave the star. A lep- 
ton fraction of 0.4 corresponds to a proton fraction of the 
order of 0.3. During deleptonization the proton fraction 
of stellar matter decreases and after total deleptonization 
at densities below 0.08 fm~'^ the proton fraction is below 
0.1. From the behavior obtained in Fig. [6] we conclude 
that during deleptonization the number of nucleons in 
the clusters decreases, the low density background gas of 
nucleons and light clusters increases and the onset of the 
clusterized phase shifts to larger densities. 

Next we analyze the properties of the clusters formed 
in the pasta phase. We consider the droplet geometry 
because of the finite size of this structure. In particular 
we discuss the effects of temperature, isospin asymmc- 
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FIG. 8: (Color online) Droplet and Wigner Seitz cell radius of 
the spherical clusters for NL3, NL3(j;p, FSUgoM and lU-FSU 



calculated at T = 5 MeV and K,, 
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try and the density dependence of the symmetry energy 
on the number of nucleons and protons in the droplets, 
on the onset density of this geometry and the transition 
density to the rod geometry. 

We denote by Adrop the number of nucleons belong- 
ing to a droplet (the type of structure that appears at 
the lowest densities in the non- homogeneous phase). We 
first check if we are including correctly the light clus- 
ters: in fact, we have to check that j4(irop > 4, to confirm 
that classifying these droplets as heavy clusters is correct. 
Secondly, we can cornpare our results to those obtained 
in other approaches, [3,13- We calculate this quantity as 

Adrop = f i?i, {p' ' p'\Yl' + y:J)] . (38) 

In a similar way, we can compute ^dropj the charge con- 
tent of a droplet. Next we show several figures including 



our results for both Adrop and ^drop for various models, 
temperatures and proton fractions. 

In Fig. [7] we display the results for different 
parametrizations: NL3 and NLBwp in Fig. [7^), NL3, 
FSUgoM and lU-FSU in Fig. W}- Here and in the fol- 
lowing analogous figures, the solid lines represent Adrop, 
the dashed ones .^drop- The onset of the droplet phase is 
characterized by a discontinuity on the number of nucle- 
ons inside the cluster: it is necessary a minimum num- 
ber of nucleons to compensate the surface energy, that is 
larger for NLSwp. In a Thomas-Fermi calculation, where 
the surface energy is calculated self-consistently we may 
expect a less discontinuous behavior. A change in the 
isovector channel of the model as in NLSwp leads to a 
large effect on the number of nucleons, increasing this 
number to more than the double. As discussed before, 
a smaller symmetry energy slope corresponds to a larger 
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FIG. 11: (Color online) Droplet and Wigner Seitz cells radii 
of the spherical clusters for lU-FSU calculated at T = 5 MeV 
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surface energy and neutrons do not drip out so easily. 
The number of nucleons obtained within NLSwp is con- 
sistent with the results of Q within a statistical model, 
based on the TMA parametrization [3a] . We should point 
out that for NL3 the separation that has been done in 
hght clusters and heavy clusters breaks down for a small 
proton fraction because, at low densities, the size of the 
heavy clusters equals the size of the light clusters. The 
models which include an isoscalar-isovcctor coupling A.„ 
present larger nuclei, the heaviest ones corresponding to 
lU-FSU. In these models, the appearance of the heavy 
clusters occurs at similar densities, larger than that ob- 
tained for NL3. This again is due to the fact that NL3 
has the largest symmetry energy slope at these densities. 
In Fig. [5] we plot the radius of the Wigner Seitz cell 
together with the droplet radius as a function of density 
for NL3, NLSwp, FSUgoM and lU-FSU. The ordering of 



the radii obtained in the different parametrizations re- 
flects perfectly the ordering of their surface tensions [cfr. 
Fig. [l]b)], that, in turn, is closely linked to the symmetry 
energy density dependence [cfr. Fig. [T]a) and Table HI] . 
NL3 has by far the smallest surface energy at Ypg = 0.2, 
while lU-FSU has the largest, see Fig. [TJd): correspond- 
ingly, NL3 has the smallest Wigner Seitz cell and droplets 
and lU-FSU the largest ones, as shown in Fig. IS] 

In Fig. [5] the ratio ^drop/^drop is plotted at a function 
of density for the different models we are comparing. We 
conclude that this ratio decreases with density and is 
model dependent. A decrease of the proton fraction of 
the clusters with density was also obtained in Q. 

The models with a smaller symmetry energy slope have 
smaller proton fractions. A smaller slope implies that 
neutrons drip out of the cluster with more difficulty giv- 
ing rise to neutron richer clusters. Also, a smaller slope 
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results in a smaller £sym.{p) above O.Tpo [— 0.1 fm~'^, 
where approximately all the curves cross, see Fig. [T]a)]. 
Since the density inside the droplets is between 0.7po 
and po J a smaller symmetry energy favors less symmetric 
clusters. 

The dependence of Adrop and Z^mp on the temperature 
for Ypg = 0.2, and on the proton fraction, for T = 5 MeV 
is plotted in Fig. [10] a) and b), respectively, for the lU- 
FSU parametrization. Decreasing the temperature in- 
creases slightly the number of nuclcons in the clusters 
and strongly decreases the onset density. The transition 
to the rod geometry seems to be temperature indepen- 
dent. 

Isospin asymmetry affects the number of nucleons in 
the cluster in a non linear way as is seen in the right panel 
of Fig. [TU] for the lU-FSU parametrization. At the onset 
of the pasta phase the size of the clusters is smaller for 
the smaller proton fractions. However, there is a faster 
increase of the cluster size with density when the pro- 
ton fraction is smaller. In particular, at the transition to 
the rod geometry we get smaller droplets in symmetric 
matter. This behavior is also obtained for FSUgoW and 
NL3a;p. However, for NL3, we get a systematic behavior, 
the smaller the proton fraction the smaller the number 
of nucleons in the cluster. In [8|, using the Skyrnie in- 
teraction SKM* with a quite low symmetry energy slope 
at saturation (L=45 MeV), a cluster size independent of 
the proton fraction was obtained. On the other hand, 
the authors of using TMA get at T = 1 MeV larger 



clusters for Y„, 



0.3 and smaller for Yr,, 



0.5, and at 



T ~ 5 MeV larger clusters for Ypg =0.3 and smaller for 
Ypg =0.1. It seems that the proton fraction is not af- 
fecting the size of the clusters in a linear way, and this 
behavior is model dependent. In fact, from Fig|6] it is 
seen that the proton fraction with the largest fraction 
of nucleons in the clusters is Ypg = 0.5. However, the 
number of nucleons in the droplets also depends on their 
size. In FiglTl]we show how the size of the Wigner Seitz 
cell and droplet radius depend on the density and pro- 
ton fraction. The number of particles in the droplets is 
strongly dependent on the these two radii. 

In Fig. [T^we show for lU-FSU how the proton frac- 
tion in the droplets depends on the temperature and on 
the global proton fraction. Temperature makes droplets 
richer in protons and isospin asymmetry reduces the pro- 
ton fraction in the droplets. In fact, temperature helps 
the evaporation of neutrons from the droplets increasing 
the proton fraction in the clusters. On the other hand, 
increasing the neutron fraction leads naturally to an in- 
crease of neutrons in the clusters. These two effects were 
also obtained in the formalism developed in Q within 
a grand-canonical ensemble approach (see left panels of 
Fig. 25 in Ref. [|). 

To complete the discussion, we now distinguish the 
contribution of each type of light cluster, i.e., we decom- 
pose Yiight- In Fig. [T2]we show the fractions of nucleons 
belonging to a's, hclions, tritons and deuterons for lU- 
FSU at T = 5 MeV and Ypg = 0.2, as a choice that 



exemplifies the results we obtain. 

Two typical sequences can be identified in the abun- 
dances of the various clusters: d, t, h, a at low baryonic 
density and d, t, a, h at the higher total baryonic density 
at which clusters are still present. Deviations from this 
behavior are found for lower temperature (T = 3.5 MeV) 
or higher proton fraction (Ypg = 0.5). In general, the 
dominant contribution comes from the deuterons. At 
very low density, the sequence reflects the ordering in 
size: the smaller clusters are more abundant. The dif- 
ference between t and h, and more precisely the finding 
that tritons are more abundant than helions, is due to 
two effects: globally there are more neutrons available 
and. moreover, the triton is more bound. This last point 
explains why we find a difference among them also at 

Ypg = 0.5. 

For a given temperature, the fraction of a-particles 
at very low density increases with density much faster 
than the other clusters: this is due to their large binding 
energy. As a consequence, the fraction of a particles may 
overcome the fraction of other type of clusters at larger 
densities. At sufficiently low temperatures, they are the 
most abundant cluster in HC matter, even for the lowest 
densities. In HPC matter, their abundance is determined 
by the density, proton fraction, and temperature of the 
background gas. 

Concerning the dissolution density of each cluster, we 
find that it increases with increasing temperature; on the 
other hand, the effect of a larger proton fraction is just 
to increase slightly the abundancy of the clusters (except 
for tritons), it docs not alter their dissolution. 

Finally, we show in Fig. [14] the average number of nu- 
cleons and the average charge of the constituents in the 
HC and HPC cases. We consider as constituents the 
protons, neutrons, a, helions, tritons, deuterons and the 
droplets. In the low density limit, the system is composed 
by free nucleons: in fact, (A) = 1 and (Z) = 0.2 in this 
limit. At intermediate densities, (A) increases slightly 
because of the formation of light clusters. At even higher 
densities, the behavior is different for the HC and HPC 
constructions. In the HC case, the clusters dissolve at 
some density and we get again (A) = 1, which means a 
gas of free nucleons. In the HPC case, droplets form, and 
the average number of nucleons in a cluster increases very 
rapidly with the density. Concerning the average charge, 
from the figure it is seen that for the droplets it also in- 
creases, but more slowly than (A): as already discussed, 
the isospin-asymmetry of the droplets increases with the 
density. 



V. CONCLUSIONS 

In the present paper we have investigated the effects 
caused by the explicit inclusion of four light clusters, 
namely, a particles, deuterons, tritons and helions, in ho- 
mogeneous and non-homogeneous nuclear matter at low 
densities. This study is particularly important for the 
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FIG. 13: (Color online) HPC case; the various contributions 
to Yiight are shown explicitly for the lU-FSU parametrization 



mctry energy. In the folio-wing we identify some of the 
trends that were found: a) when temperature increases, 
the pasta phase appears at higher densities, an effect al- 
ready seen without the inclusion of the light clusters; b) 
the pasta phase low density gas, including light clusters, 
has a larger number of particles at larger temperatures; 
c) larger abundances of light clusters occur at densities 
below the pasta phase onset and their fractions are larger 
for more symmetric matter. Before the onset of the pasta 
phase their abundances decrease with temperature, and 
the contrary occurs in the pasta phase; d) within the 
coexisting-phases method adopted in the present work, 
the onset of the pasta phase occurs at lower densities for 
more asymmetric matter. We believe this is due to the 
zero thickness surface approximation and a lack of self- 
consistency in the calculation of the surface properties. 
An opposite trend was obtained within a Thomas Fermi 
approach 3J|. However, if we also consider the light clus- 
ters, the onset of clusterization starts at lower densities 
for the larger proton fractions. 
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understanding of the composition of the inner crust of 
protoneutron stars. 

We have chosen to calculate the above mentioned ef- 
fects with four different parametrizations of the NLWM 
because they are characterized by different symmetry en- 
ergy density dependencies, as seen from the values of 
their slopes L in Table II. 

We have checked that the influence of the light clusters 
in the free energy per particle is only noticeable at very 
small densities (up to 0.0025 fm~'^) both in homogeneous 
and pasta phase matter. The inclusion of light clusters 
lowers the free energy and their effect is smaller in the 
pasta phase range. 

We have analyzed the fractions of nucleons at different 
temperatures and different proton fractions. The results 
are model dependent, as expected, and some of the dif- 
ferences are related to density dependence of the sym- 



The number of nucleons and the number of charged 
particles inside the pasta droplets were calculated. It was 
shown that: a) models with a smaller symmetry energy 
slope have larger clusters, inside a larger Wigner-Seitz 
cell, with a larger number of particles and smaller proton 
fraction; b) the number of nucleons in the droplets has 
no linear relation with the global isospin asymmetry of 
matter and density, showing a strong dependence on the 
properties of the symmetry energy; c) the number of nu- 
cleons in the clusters decreases if temperature increases; 
d) the fraction of protons in the clusters decreases with 
density, decreases for more asymmetric matter and in- 
creases when temperature increases. 

We have shown that the composition and structure of 
the pasta phase is quite sensitive to the symmetry en- 
ergy behavior at low densities. It would be interesting 
to study the transport properties and neutrino opacity 
in the pasta phase with light clusters, because these are 
important quantities in the cooling mechanism of pro- 
toneutron stars. 
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Appendix 

The surface tension is obtained by fitting the Thomas- 
Fermi results with the formula 



cr(x, T) = a{x) [1 - a{T)xT - b{T)T^ - c{T)x^T] 
where 



(39) 



a{x) = (Toexp (— cria;'^")(l + flia; + a2a; 



a^x + a/ix + a^x + aga; 

2 I „ ^3 



a{T) = flo + aiT + a2T + a^T 
h{T) = aa + aiT + a2T^ + a^T^ 



c{T) = gq + aiT + a2T + asT' + a4T 






(40) 
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Clearly, a{x) is the surface tension at T = 0, and do is 
its value at T = for symmetric matter. Notice moreover 
that there are a few more terms in these parametrizations 
as compared to those used in Ref. |19| . 
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TABLE III: Surface tension parameters fitted within the 
Thomas-Fermi approximation for NL3, NL3a;p, FSUgoM and 
lU-FSU parametrizations. The coefficients are for T in MeV 
and (To is in MeV/fm~^. 
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NL3 


a{x) 


a{T) 


b{T) 


c{T) 


O"0 


1.12307 


- 


- 


- 


0-1 


20.7779 


- 


- 


- 


Oo 


- 


0.0121222 


0.00792168 


- 


ax 


-5.84915 


0.01664 


-8.2504x10"^^ 


- 


a2 


138.839 


-0.00137266 


-4.59336x10"'' 


- 


as 


-1631.42 


4.0257x10"^ 


-2.81679x10"'^ 


- 


a4 


8900.34 


- 


- 


- 


as 


-21592.3 


- 


- 


- 


ae 


20858.6 


- 


- 


- 


NL3a;p 


a{x) 


a{T) 


b{T) 


c{T) 


cro 


1.12013 


- 


- 


- 


cri 


14.0774 


- 


- 


- 


Oo 


- 


-5.80451x10"^ 


0.00725961 


-0.00259094 


Ol 


-2.15376 


0.0233833 


0.000318409 


-0.053756 


02 


57.8455 


-0.00507732 


-0.000104941 


0.0114598 


as 


-431.365 


0.000490863 


1.19645x10"^ 


-0.000354375 


a4 


1854.81 


-1.59473x10"^ 


-7.19099x10"'^ 


-4.76451x10"^ 


as 


-3653.96 


-7.55062x10"** 


1.62087 xlO~* 


2.28389x10"^^ 


ae 


3214.82 


- 


- 


- 


FSUgom 


a{x) 


a{T) 


6(T) 


c{T) 


CTO 


1.1223 


- 


- 


- 


o-i 


-1.45717 


- 


- 


- 


OO 


- 


-0.0133789 


0.00773356 


0.0408077 


Ol 


-3.17729 


0.0330912 


-0.000240406 


-0.0971609 


02 


-9.5121 


-0.00786564 


4.52523x10"^ 


0.0195288 


OS 


70.5609 


0.000902286 


-7.64893x10"'' 


-0.00140166 


04 


-155.641 


-4.84828x10"^ 


5.33346x10"^ 


4.97386x10"^ 


OB 


154.691 


9.56728x10"^^ 


-1.45394x10"* 


-1.20803x10"*^ 


Og 


-58.9476 


- 


- 


- 


lU-FSU 


a{x) 


a{T) 


6(T) 


c{T) 


O"0 


1.16473 


- 


- 


- 


CTl 


-0.659167 


- 


- 


- 


OO 


- 


0.00404325 


0.00767923 


0.0066774 


Ol 


-2.25482 


0.00828207 


-8.58068x10"^ 


-0.0514285 


02 


-5.64237 


-0.00153301 


4.43918 xlO"^ 


0.00949505 


OS 


37.8471 


7.26763x10"^ 


-5.44453x10"'^ 


-0.000427613 


04 


-81.6617 


- 


- 


- 


OB 


81.2696 


- 


- 


- 


ae 


-31.0227 


- 


- 


- 



